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Abstract 

Semiconductor spintronics will need to control spin injection phenomena in the non-linear regime. 
In order to study these effects we have performed spin injection measurements from a dilute mag- 
netic semiconductor [(Zn,Be,Mn)Se] into nonmagnetic (Zn,Be)Se at elevated bias. When the ap- 
plied voltage is increased to a few mV we find a strong decrease of the spin injection efficiency. 
The observed behavior is modelled by extending the charge-imbalance model for spin injection to 
include band bending and charge accumulation at the interface of the two compounds. We find that 
the observed effects can be attributed to repopulation of the minority spin level in the magnetic 
semiconductor. 

PACS numbers: 72.25.Dc, 72.25.Hg, 81.05.Dz 



Much of our understanding of the physics of electrical spin injection from a ferromag- 
net into a normal metal derives from a current-imbalance model proposed by van Son0 
and, independently, by Johnson and Silsbee0. In this model, the current conversion is 
driven by a splitting of the Fermi levels for majority- and minority-spin electrons at the 
ferromagnet-normal metal interface. The concomitant reduction of the conductivity of the 
normal metal due to the suppression of a spin-channel leads to a discontinuity of the aver- 
age electrochemical potential, which is sometimes referred to as the spin-induced boundary 
resistance. The model was put on a solid (Boltzmann) basis by Valet and Fert|§ and proved 
very useful for describing important spin transport phenomena including CPP-GMR (giant 
magnetoresistance in the current-perpendicular-to-plane geometry) Q . 

In a further application, we used the model to describe spin injection into 
semi conductors 0, revealing the importance of a conductance mismatch between ferromag- 
netic metals and semiconductor. The conductance mismatch keeps the splitting of the 
spin-dependent Fermi levels at minimum, thus precluding spin injection. In the meantime, 
several methods of avoiding the conductivity mismatch have been proposed. However, by 
far the most robust route towards spin injection to date|| |7] is the use of dilute II- VI mag- 
netic semiconductors (DMSs) that exhibit the giant Zeeman effect 0, have a conductivity 
comparable to that of non-magnetic semiconductors, and can boast spin polarizations close 
to 100 % at a small applied magnetic field. 

Recently, we used a simple DMS - non-magnetic semiconductor (NMS) heterostruc- 
ture (consisting of (Zn,Mn,Be)Se as DMS and lattice matched (Zn,Be)Se as NMS) to 
demonstrate the field dependence of the spin- induced boundary resistance, which increases 
with field as the magnetization of the paramagnetic DMS increases. The experiments re- 
ported in Ref. |7j were all done in the regime of linear response, where the current-imbalance 
model is appropriate. 

However, spin injection experiments in semiconductors allow one to very easily enter 
the regime of non-linear response, where corrections to the current-imbalance model are 
necessary. These corrections will be manifest in any future spintronic semiconductor device, 
and it is thus of importance to identify them correctly. 

In this paper we report on spin-injection measurements in the non-linear regime. More 
specifically, we experimentally demonstrate the strong effects that band-bending and charge- 
accumulation have on the non-linear transport. We present a modelling of the observed 
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phenomena which generalizes the charge-accumulation model to include the limited screening 
available in semiconductor structures, and shows good agreement with the experiments. It 
should be noted that the effects discussed here involve non-linearities influencing the spin 
polarization at the DMS/NMS interface, and thus are of a fundamentally different nature 
than the drift-induced effects discussed recently by Yu and Flatte0 (which occur at higher 
electric fields). 

The devices on which the experiments were carried out consisted of an all II- VI semi- 
conductor heterostructure fabricated by molecular beam epitaxy. The heterostructure con- 
sisted of three semiconductor layers. From bottom to top these layers are a non mag- 
netic n-type Zno.97Beo.03Se layer (thickness 500 nm, n ~ 10 19 cm -3 ), a dilute magnetic 
Zno.89Beo.05Mno.06Se layer (thickness 100 nm, n « 5 x 10 18 cm -3 ) which acts as a spin 
aligner, and a top layer of 30 nm highly n-doped ZnSe (n ~ 2 x 10 19 cm -3 ). The latter 
was grown to ensure good quality ohmic contacts and was covered in-situ with aluminium 
without breaking the UHV. 

In the Al-layer, 200 x 250 /mi contact pads were defined by optical lithography and wet 
chemical etching. These pads were used as a mask for a subsequent wet etching step through 
the ZnSe and the DMS down to the Zno.97Beo.03Se. A mesa was then defined by etching 
down to the substrate, leaving only two contacts and the transport layer in between. The 
resulting sample is schematically shown as the inset in Fig. [I]. 

The samples were immersed in a magnetocryostat and their transport properties were 
determined at temperatures of 1.6 K, 3 K, 4.2 K and 6 K. The magnetoresistance of the 
devices was measured using dc techniques and a quasi four probe geometry, in which the 
wiring resistance of the setup was excluded, while the contact resistance of the device was still 
part of the measured resistance. [|10| For bias voltages Vbi as in the regime of linear response 
(300/iV or less) the device showed a positive magnetoresistance. Fig. [I] plots the relative 
magnetoresistance AR/R for a sample with a distance xq = 10 /xm between the contact 
pads, taken at 1.6 K, where the zero-field resistance R = 420 Q. As described in Ref. J7|, 
the magnetoresistance results from the increase of the spin-induced boundary resistance 
with magnetic field at the two DMS/NMS interfaces in the samples. All data discussed here 
were taken on the same sample as that for which the low-temperature magnetoresistance is 
shown in Fig. [I]; we verified that the effects discussed occur in samples with varying doping 
concentrations and dimensions. We found experimentally that both R and the saturated 
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magnetoresistance AR/R w 0.25 are independent of temperature in the investigated range. 

The main experimental result of this paper is shown in Fig.|2|(a) . When the applied voltage 
is increased, a pronounced and very rapid drop of the magnetoresistance is observed, reducing 
the effect by two or more orders of magnitude on applying a voltage of around 10 mV across 
the junction. The experimental data in Fig. 2a were taken starting from three different 
values of AR/R (i.e., at different values of the magnetic field B) in the linear response 
regime (i.e., AR/R « 0.05, 0.1 and 0.15, respectively), at the four different temperatures 
mentioned above. Obviously, the non-linearities show a marked temperature dependence. 
Moreover, the horizontal axis displays the voltage drop over the junction Vj, rather than 
the total bias V^as (roughly, Vj ~ 0.15Vbi as ), for reasons of comparison with the modelling 
discussed below, and because it more clearly illustrates the energy scales involved in the 
non-linearities. We will detail below how Vj is defined and can be calculated from Vbias- 

The drop of the magnetoresistance can be understood if we combine the model for diffu- 
sive spin polarized transport with the band structure of the semiconductor heterostructure. 
When a current is driven from a spin polarized material into a non polarized material, the 
electrochemical potentials for spin up (/r) and spin down (/x^) split at the interface. In 
linear response, JTT| the length scale of this splitting is given by the spin scattering length of 
each material. The situation is depicted in Fig. 3, where the Zeeman-split conduction band 
(full drawn lines) and relevant potentials (dash-dotted lines) of the DMS are shown of the 
left side of the figure. The interface is indicated by the dotted vertical line at x — 0, with 
the NMS in the right half of the plane. The splitting of (jj,*) and (/i^-) is the driving force 
which leads to a spin polarized current in the non magnetic material. Because the conduc- 
tivities for spin up and spin down are equal in the NMS, only a difference in the derivative 
of the electrochemical potential can lead to different currents in both spin channels. Since 
the electrical potential must be equal for both spin directions, this difference can only be 
introduced through the chemical potential, i.e., by spin accumulation. Spin injection thus 
leads to a potential drop at the interface which drives the spin conversion. This voltage 
drop, which may alternatively be regarded as a spin-induced boundary resistance (since it 
represents an extra voltage drop in the sample in absence of spin injection), is indicated in 
Fig. 3 by the potential difference at x = between the thin drawn lines, denoted /i^ MS and 
A*nms> that extrapolate the voltage drop in the bulk of the material towards the interface. 

While in the NMS the splitting of the Fermi levels is symmetrical because the conduc- 
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tivities for spin up and spin down electrons are identical, in the DMS, the splitting for 

the majority- (c T = ^dms(°) ~ ^dms( )) and the minority- (c l = ^Sms(°) - /4)ms( )) s P in 
electrons can, in a one-dimensional model, be expressed as[]7|, [12] 



C T c l = A NMS Iftjft ± 1) ,y 

0~NMS (I + e _X NMS + 2^^e _ *NMs') + ^NMS£DMS (I _ QZ) 
V XO I "NMS^DMS V r ' 

where in the numerator the plus(minus) sign applies to c 1 (c*), respectively. In Eq. 1, 
Adms, Anms, o"dms, 0"nms, are the spin flip length and the conductivity in the DMS and the 
NMS respectively, x is the spacing between the contacts, I is the current and ft is the 
degree of spin polarization in the bulk of the contacts. Note that and are defined 
setting /iD MS (0) as reference level for the energy scale, i.e. /iD MS (0) = 0. For the potential 
drop AU at the interface, and the resulting magnetoresistance we simply have: 



Atf = K) V ms(0) - /4 V MS (0) = (c T + ci)/2, (2a) 

AR = AU/I. (2b) 

Eqs. 1 and 2 are quite general and describe spin injection in metals as well as semicon- 
ductors - but only in the linear regime. The magnitude of the Fermi-level splitting (and 
thus of AU) is different for different types of junctions: since the spin-polarized current is 
driven solely by the spin accumulation, the Fermi-level splitting has to be of the order of 
the current imbalance between the spin channels times the resistivity of the normal metal. 
When the magnet and the non magnet are both metals, as, e.g., in CPP GMRQ, the split- 
ting remains small even for high currents due to the high conductivity. When the magnet is 
metallic and the non-magnet is a semiconductor, the splitting is determined by the material 
with the highest conductivity, and thus is also small. This effect leads to negligible levels of 
spin injection ||, and is also known as the conductance mismatch. When magnet and non 
magnet are both semiconductors as in the large magnetoresistance effect discussed here, 
the splitting can easily be in the range of mV. 

In a magnetic metal the Fermi energy is typically of the order of a few eV. A Fermi-level 
splitting of a few /xV at the interface will thus not influence the properties of the magnet, 
and one expect the linear regime to hold up to quite high current densities. In a DMS the 
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situation is completely different. The Fermi energy, the Zeeman splitting, and the Fermi 
level splitting are in all the range of mV, and non-linear effects are to be expected when bias 
voltages of similar magnitude are applied. 

This situation, applied to the present experiment, is pictured in some detail in Fig. |3|. 
The conduction band of the NMS is some tens of mV below the conduction band of the DMS. 
The conduction band level of the DMS is split by the Zeeman splitting into two subbands, 
E™ s and E^ s . From previous spin injection experiments^, [7|] and from spin flip Raman 
scattering we know that the DMS is fully spin polarized at low temperatures and moderate 
magnetic fields, which recent band structure calculations understand as resulting from the 



formation of an impurity band. [ 13 1 This implies that the Fermi energy is situated above the 
lower and at least a few mV below the upper Zeeman level. 

As discussed above, spin injection will lead to the occurrence of a 'built-in potential' AU 
at the interface. This is an actual electrochemical potential step (i.e., not spin-dependent). 
In order to preserve both charge conservation and the band offset at the junction, AU has to 
be compensated by band bending and charging at the interface. In Fig. 3, this is indicated 
by the dashed lines emanating from E^f 13 and E^ s . (In principle, one also expects band 
bending at the NMS side of the junction. For clarity, we have not included this in Fig. 3, 
nor in the modelling we describe below. Its inclusion is straightforward.) 

Thus modified, Fig. 3 satisfies all boundary conditions imposed by the band structure 
and the continuity equations for the spin-polarized currents at the interface and for the 
total current throughout the device. However, it is now obvious that the properties of 
the DMS can be seriously affected by the band bending. At the interface, the majority spin 
electrochemical potential ^ now approaches the upper Zeeman level -E^?| MS . From the strong 
temperature dependence observed for the LMR effect in Ref. we know that this energy 
gap is small enough to strongly influence the population of the upper Zeeman level and thus 
the bulk spin polarization (5 in the DMS, even at low temperatures. (3 being close to one, 
however, is a prime prerequisite for injecting a highly spin polarized current into the NMS. 
We can thus expect the spin polarization injection (and thus the LMR effect) to collapse as 
soon as the band bending starts to reduce (3. 

Since j3 and AU depend on each other, a modelling of the phenomena as a function of 
applied bias voltage should be done in a self-consistent manner. However, one can formulate 
a mathematically consistent description when starting from a given value of AR/R in the 
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linear response regime (this is the main reason for presenting the data emanating from the 
same AR/R value in Fig. 2.) We first use Eqs. 1 and 2 to calculate the bulk polarization j3 
in the DMS. In the linear regime, the bulk value of (3 equals /3(x = 0), the spin polarization 
at the interface. Assuming Boltzmann statistics, we then directly have the energy splitting 
between E%f s (0) and /i T (0) from 

p{x = 0) = tanh((£D MS (0) - AU - rf(0))/2k B T), (3) 

where of course AU = for infinitesimally small bias. For modelling the dependence 
on Vu&s, we gradually increase AU (note that at this point we assume all band-bending to 
occur in the DMS), calculate the reduced /3(x = 0) using Eq. 3, and substitute this value 
for the bulk polarization in Eqs. 1 and 2 to calculate AR/R. At the same time, AU can be 
converted in a voltage drop across the junction, Vj. 

The latter quantity is conveniently accessible for comparison with the experiment. This is 
because AR — > 2Anms/ ct nms f° r B, (x /Anms) oo, as can easily be verified from Eqs. 1 and 
2. Experimentally, we have (within our 1-dimensional modelling) onms — 2.5 x 10 _4 fi~ 1 cm, 
yielding Anms — 1.25/mi. For both the experimental (Fig. 2a) and theoretical (Fig. 2b) 
plots of the non-linear behavior, we may now calibrate the voltage axis according to Vj = 

IAR + Vbias(A N MsM))- 

As to the remaining parameters in our model, we have from experiments on single DMS 
layers that 0dms> converted to the 1-dimensional model equals 1.0 x 10 f2 _1 cm. The only 
free parameter now left in the model is Adms- Since there is no easy method to measure 
Adms; and moreover its magnetic field dependence is unknown, we simply have set the ratio 
Anms ^dms/^nms Adms in Eqs. 1 equal to 1, yielding A DM s = 0.5/im. 

The modelling of the band bending effect as described above leads to the plots shown in 
Fig. 2b. We find that indeed a few mV of voltage drop across the junction are enough to 
reduce the spin polarization of the injected current to almost zero. The computed curves 
closely resemble the experimental results in shape, magnitude, voltage range, and tempera- 
ture dependence, and confirm the validity of the modelling discussed above. 

At this point, we should address the drift effects introduced by Yu and Flatted , that 
also can induce a reduction of AR/R in our experiments. For the highly (i.e., above the 
metal-insulator transition) doped samples used here, it is easy to show that drift effects can 
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only be expected for much higher electric fields than those needed to actually observe the 
non-linearities. Moreover, within the drift model one would not expect any temperature 
dependence for degenerate semiconductors, again in contradiction with the experiments. 
We do, however, have some indications that drift effects actually occur in certain ranges of 
operation. Note, e.g., that AR/R is not completely suppressed for high Vj in the modelling 
curves of Fig. 2b, while experimentally AR/R vanishes entirely. The drift mechanism can 
actually explain this behavior, be it for fields that are an order of magnitude larger than 
those needed to induce the initial non-linearities. 

In conclusion, we have shown that when spin injection into semiconductors is used beyond 
the regime of linear response, the splitting of the electrochemical potentials can influence the 
band bending in spin-injecting junction. For finite 'built-in potentials' AU it is important 
that the energy difference (E^f 13 — /^(0)) is kept sufficiently large to ensure large spin 
polarization at the interface. This observation may seriously limit the operational voltage 
of spintronic devices and forbid applications in power electronics. Appropriate tailoring of 
the band structure may be possible to circumvent the problems described here. 
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fig 1 Schmidt et al. 

FIG. 1: Insert: Spin injection device used in the experiment consisting of a non-magnetic semi- 
conductor layer with two DMS top-contacts. The graph gives the resistance change AR/R versus 
magnetic field B. 
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fig2b Schmidt et al . 

FIG. 2: (a) Experimental and (b) theoretical non-linear magnetoresistance AR/R data plotted 
over voltage drop Vj. To facilitate comparison between experiment and theory, curves are plotted 
starting at several fixed values of AR/R (obtained by carefully adjusting B), for temperatures of 
1.6, 3, 4.2 and 6 K. The parameters involved in the modelling are discussed in the text. 
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fig3 Schmidt et al . 

FIG. 3: Diagram of the band bending at the spin- injecting DMS/NMS interface. A£" c t denotes 
the location of the conduction band offset between E^^ s and Eq MS when bend-bending is taken 
in to account; all other symbols are discussed in the text. Note the discontinuity between /uf^ MS 
and /uj^g at the junction (x = 0), which is the potential difference AU in Eqs. 2. 
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